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We introduce a set of generalized slave-particle models for extended Hubbard models that treat lo¬ 
calized electronic correlations using slave-boson decompositions. Our models automatically include 
two slave-particle methods of recent interest, the slave-rotor and slave-spin methods, as well as a 
ladder of new intermediate models where one can choose which of the electronic degrees of freedom 
(e.g., spin or orbital labels) are treated as correlated degrees of freedom by the slave bosons. In 
addition, our method removes the aberrant behavior of the slave-rotor model, where it systemati¬ 
cally overestimates the importance of electronic correlation effects for weak interaction strength, by 
removing the contribution of unphysical states from the bosonic Hilbert space. The flexibility of 
our formalism permits one to separate and isolate the effect of correlations on the key degrees of 
freedom. 


I. INTRODUCTION 

One of the long-standing areas of interest in condensed 
matter physics, particularly that of complex oxides, is 
that of the Mott metal-insulator transition 1 . Generically, 
within a Hubbard model framework, as the strength of 
localized electronic repulsions is increased, the electrons 
prefer to be localized on atomic sites and inter-site hop¬ 
ping is suppressed, and at a critical interaction strength 
the system becomes an insulator. An illustration of the 
rich behavior that can occur in such systems is the Or¬ 
bital Selective Mott Transition (OSMT) whereby only a 
subset of localized orbitals become insulating (localized) 
while the remainder have metallic (extended) bands. An 
example is provided by quasi-two-dimensional Mott tran¬ 
sition in the Ca 2 - x Sr x Ru 04 family, where the Mott 
metal-insulator transition and its magnetic properties^ 
at the critical doping x = 0.5 show a coexistence be¬ 
tween a magnetic susceptibility that shows a Curie form 
for S' = 1/2 and a metallic state. Anisimov et aZP have 
used DFT+DMFT to explain this situation in terms of 
an OSMT in which one Ru 4d orbital is localized, while 
the other continues to present metallic behavior. 

The present day workhorse for ab initio materials mod¬ 
eling and prediction, Density Functional Theory (DFT), 
is fundamentally based on band theory and is unable 
to describe such transitions (without symmetry break¬ 
ing of the electronic degrees of freedom: e.g., spin or 
orbital polarization). To this end, Hubbard model based 
methods such as Dynamical Mean Field Theory (DMFT) 
and DFT+DMFl®^ have been developed to include lo¬ 
calized correlation effects in electronic structure calcu¬ 
lations. However, DMFT-based methods are computa¬ 
tionally expensive and typical present day calculations on 
real materials are generally restricted to treating a few 
correlated sites. Therefore, it is of significant interest to 
have computationally inexpensive, but necessarily more 
approximate, methods that include correlations and can 


permit one to rapidly explore the qualitative effects of 
electronic correlations. 

One set of such approximate methods that have been 
of recent interest are slave-particle methods. Slave-boson 
methods have a long background in condensed matter 
theory for analytical treatments of correlations typically 
in the limiting case of infinite correlation strengtfPtS. 
Kotliar and Ruckenstein 11 lused a slave-particle represen¬ 
tation to treat Hubbard-like models at finite interaction 
strength, which found applications in the realm of high- 
temperature superconductivity 1 ®. Further, Kotliar and 
Ruckenstein’s model has been generalized to multi-band 
modeli® 16 where, e.g., the effects of multiple orbitals, 
orbital degen eracy, and the Hund’s interaction have been 
studiedi^ES However, the approach of Kotliar and Ruck- 
enstein, and its various extensions, require a large num¬ 
ber of bosonic slave particles: one needs one boson per 
possible electronic configuration on a correlated site. 

For this reason, more economical slave-boson repre¬ 
sent at ionshave been of significant interest. Florens and 
GeorgesP 7 used a single “rotor” slave-boson per site 
that describes the total electron count on each site in 
a computationally economical manner. The slave-rotor 
method has successfully predicted a number of electronic 
phases of nickelate heterostructures 19 which was a dis¬ 
tinct improvement over previous studies. However, a 
rotor-like description is not orbitally selective as it can 
only describe the total electron count on a site and not 
its partitioning among inequivalent orbitals on that site. 
An alternative slave-particle approach is to treat each 
localized electronic state (i.e., a unique combination of 
spin and orbital indices) with a slave boson: this “slave- 
spin” approach automatically handles orbital symmetry 
breaking and can predict OSMT^212U. Recently, it has 
been applied to predict key physical characteristics in 
iron superconductors 2 ®. 

In this work, we introduce a generalized framework 
for slave-particle descriptions. This produces a ladder 
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of correlated models, and the slave-rotor and slave-spin 
are automatically included as two specific cases. Our 
approach does not require any physical analogies to cre¬ 
ate the slave bosons (e.g., a quantum rotor or angular 
variable to motivate the slave-rotor or a pseudo-spin to 
motive the slave-spin) and works directly in the occu¬ 
pation number representation. In our approach, one can 
choose which degrees of freedom are treated as correlated 
degrees of freedom (e.g., total electron count on a site, 
electron counts in each orbital, electron count in each spin 
channel, etc.) so that we can isolate the effect of correla¬ 
tions on the separate degrees of freedom in a systematic 
manner. Section [TT] presents our general formalism, how 
it builds upon previous models, as well as gives a few 
examples of models that can be built within this frame¬ 
work. Section III is devoted to tests of possible models 
built within this formalism in a mean-field approach at 
half-filling within a one-band and a two-band model in 
order to compare our results with those of previous work 
as well as to better understand the role of the different 
terms in an extended Hubbard model within our formal¬ 
ism. In Section |IV| we conclude this paper and discuss 
possible new avenues for researchers to use this method 
and possible developments of it in predicting properties 
of correlated materials. 


II. THE GENERALIZED SLAVE-PARTICLE 
REPRESENTATION 


In this section we introduce our generalized slave- 
particle representation. In appropriate limits, our ap¬ 
proach reproduces previous frameworks such as the slave- 
rotor and slave-spin methods. One utility of our ap¬ 
proach is that it allows us to unite these two, as well 
as other intermediate models, into a single slave-particle 
methodology. A variety of slaves-particle models can be 
investigated and compared so that one can isolate which 
specific correlated degrees of freedom are critical for de¬ 
scribing a specific physical problem. 


A. Extended Hubbard model 

The general correlated-electron Hamiltonian we con¬ 
sider is an extended Hubbard model given by 


Umi'm'a is the spin-conserving hopping element connect¬ 
ing orbital imcr to i'm!a. The d are canonical fermion 
annihilation operators. We take the interaction term to 
have the standard Slater-Kanamori fornPH 
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The first and second term stem from Coulombic repulsion 
terms between the same spatial orbital (U) and different 
spatial orbitals ( U'). The third term is Hund’s exchange 
between different orbitals of the same spin with strength 
J. The fourth term contains the intrasite “spin flip” and 
“pair hopping” terms. The index d is the spin opposite 
to a. The subscripts i on the U, U' and J parameters 
denote the fact that each correlated site can have its own 
set of parameters; however, to keep indices to a mini¬ 
mum below, we suppress this index. The various number 
operators are 

ima 
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For what follows, we keep in mind that due to the fact 
that nj ma = hi mai the Hund’s term in Hl nt can be rewrit¬ 
ten in an equivalent form to give 
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The index i ranges over the localized sites of the sys¬ 
tem (usually atomic sites), m ranges over the localized 
spatial orbitals on each site, a denotes spin, H\ nt is the 
local Coulombic interaction for site i detailed further be¬ 
low, tima is the onsite energy of the orbital ima , and 


B. Spinons and slave bosons 

The interacting Hubbard hamiltonian is impossible to 
solve exactly and even difficult to solve approximately. 
Part of the difficulty comes from the fact that we have 
interacting fermions which have both charge and spin 
degrees of freedom. Following well-known ideas in slave- 
boson approachetP®^, one separates at each site the 
fermionic degrees of freedom from the charge degrees of 
freedom by introducing a bosonic “slave” particle on that 
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site. The boson is spinless and charged, and one also has 
a remaining neutral fermion with spin termed a spinon. 

With spinons denoted by / operators and slave bosons 
by O operators, we define 

dimer fimaOia ( 3 ) 


Equations ([3]) and Q mean that the number of spinon 
and slave particles on each site must track each other 
because they are annihilated and created at t he sa me 
time. Hence, following ideas from prior work, 1 ' 20 l we 
enforce constraints to ensure the particle numbers track 
each other. We enforce the constraint 
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Requiring the / to be fermionic field operators in turn 
requires the O operators to obey bosonic commutation 
relations. We note that while the / spinon operators are 
standard fermionic Fock field operators, the bosonic O 
operators are generic and ad hoc: there is no assumption 
or requirement that the O be bosonic operators for a Fock 
space (and in general they are not of that variety). 

The index a is part of our generalized notation that 
permits us to unify many slave-particle models. The 
meaning of a depends on the type of model chosen, as 
we will show in detail below with a variety of examples. 
The index a refers to a subset of the complete set of mo 
indices that belong to a site i. For example, if we use 
an 0(2) slave-rotor model for the correlated orbitals on 
an sit e! 17 * 18 * where O = e~ lS and 8 is the phase angle of 
the 0(2) rotor, then a is nil: Oi a = Oi. Namely, we 
have a single slave particle on each site i that tracks the 
total number of particles on that site. At the opposite 
limit, we can have a unique slave boson for each mo (the 
“slave-spin” methocP^lSlJ) ; so that a = mo. 

We work directly in the number representation and 
introduce a number operator for the slave particles (this 
is a generalization of the angular momentum operator for 
the slave-rotor approach or the S quasi-spin of the 
slave-spin representation). The minimum and maximum 
allowed particle numbers are N m in and N max so that 

Ni a = diag (A min > -A m in T 1 , • • • , -A , m ax 1 5 A max ) • (5) 


This operator simply keeps track of the number of slave- 
particles in each slave mode ia , and the minimum and 
maximum allowed occupancies depends the slave model 
we choose as discussed below. 

Since we have introduced new degrees of freedom and 
enlarged the Hilbert space of the problem, it is neces¬ 
sary to enforce constraints so that one avoids consid¬ 
ering “unphysical states” that have no correspondence 
to those in the original problem. The original Hilbert 
space is spanned by kets of the form |{nimo-}) in the 
occupancy basis of the di ma operators. The enlarged 
Hilbert space is spanned by kets of the product form 
\{ni m a}}f\{N ia }} s where the subscripts label spinon and 
slave sectors. Within this enlarged space, there is a 
subset of “physical states” that correspond to the orig¬ 
inal kets. The first part of the correspondence is make 
the fermionic occupancies {ni mcr } of the original electron 
counts (d) and the spinon counts (/) identical. Hence, 
the real question is which Ni a are physically allowed. 


£ flmJimM = N ia \4>) ( 6 ) 
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on the kets \(f >) in the enlarged space. Only these kets are 
physically allowed in the exact description of the system, 
and they span the physical subspace. We note that the 
constraint is on the allowed kets and not on the oper¬ 
ators: the operators act in the extended Hilbert space 
that includes physical and unphysical unphysical states, 
and the f or O operators acting alone can move us from 
a physical state to an unphysical state. 

Enforcing the constraints of Eq. © ensures that only 
physical states that are in one-to-one correspondence 
to the original states are considered in the extended 
spinon+slave boson Hilbert space. (This is the same 
idea as Ref. (H3 where the 0(2) rotor angular momentum 
operator L has been replaced by N ia .) When the con¬ 
straints are obeyed, only physically allowed occupancies 
for the bosonic operators are relevant which is the same 
as setting N m i n = 0 and N max to the maximum allowed 
occupancy in the slave sector on the site. For example, 
for a system of 3 spatial orbitals on a state (i.e, three 
choices of m), N max = 6 if we only wish to count total 
numbers of electrons using the slave bosons in which case 
a is nil. Or, if we want to count up and down spin elec¬ 
trons separately for this site, then a is the same as o and 
Nmax = 3 for each choice of a. However, as explained 
below, when one does approximate calculations, choosing 
N m i n and N max that differ from these values allows for 
the creation of different types of models and offers some 
technical advantages. 

For completeness and clarity, Appendix [A] provides ex¬ 
plicit examples of the enlarged Hilbert spaces and various 
choices of a and explains how the action of the di mry and 
the fimaOia operators are identical on the physical sub¬ 
space of states. 


C. Slave operators and Hamiltonian 


To reproduce the standard behavior of the annihilation 
operator where only ph ysical states are allowed in the 
spectrum of the operator^fillil (also see the Appendix [A|, 
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However, if nima = 0, then the action of fi ma will destroy 
the total state | <j>) regardless of what O a may do, so for 
this case we have an undetermined situation: 

Oia\Ni a = 0) s = undetermined . (10) 

Following the same logic for the creation operators yields 

dl\N ia ) s = \N ia + l) s (11) 

until we reach the ceiling N^ a = N max when we have a 
similar indeterminacy 

0\ a \N max ) s = undetermined. 

Putting this all together, the slave boson operator O ia 
in the number basis must have the form 


D. Decoupling spinons and slaves 

Up to this point, our considerations have been for an 
exact solution of the interacting Hamiltonian which is im¬ 
possible in practice. To make progress, in slave-particle 
approaches one splits the problem into two separate and 
simpler pieces that are connected to each other via self- 
consistent averages of the relevant operators. Specifically, 
the ground-state wave function of the system is approx¬ 
imated by a simple, separable product state |’3//)| < F S ) 
where |'P_f) is a spinon wave function and |$ s ) is a 
slave boson wave function. Since the spinons and slaves 
are now decoupled and their number fluctuations are no 
longer locked in step, the constraint of Eq. © can only 
be satisfied on average, 
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where Ci a is at this point an undetermined constant that 
we are free to choose. Below, we will use this freedom to 
ensure that we reproduce a desired non-interacting band 
structure at zero interaction strength (when H- nt = 0). 

Substituting the spinon and slave operators into the 
original extended Hubbard Hamiltonian gives the follow¬ 
ing form, which for the moment we specialize to the sym¬ 
metric U' = U,J = 0 case to keep the logic simple (the 
more general cases are enumerated further below): 
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For the onsite tima terms, we have replaced 
flmafimaO\ a d ia by the simpler because 

even though d\ a O ia is not necessarily identity (unless 
Cia = 1), the two set of operators act identically on all 
the physical states of interest due to the fact that fi ma 
annihilates the state with zero particles. 

The point of introducing the slave degree of freedom 
is that they track the number of electrons in the vari¬ 
ous orbitals and spin states on each site and the inter¬ 
action Hamiltonian of Eq. ([2| is essentially determined 
by these numbers. Hence, one can write the most impor¬ 
tant parts of the interaction Hamiltonian solely in terms 
of the bosonic slave operators. 


where the / and s subscripts denote averaging over the 
spinon | f) and slave boson |$ s ) ground state wave func¬ 
tions, respectively. 

A priori, one can make only few statements about 
when such a decoupling scheme is expected to be a 
good approximation. First, at zero interaction strength 
(U = J = 0), the decoupled approach can reproduce the 
non-interacting band structure since the spinons alone 
can do this. Second, if number fluctuations about the 
averages are small so that imposing Eq. (13) is close 
to imposing Eq. © , then we expect the decoupling to 
work well in describing ground-state averages; examples 
of phases with small or zero number fluctuations are nar¬ 
row or Mott-like insulating bands. Third, in infinite di¬ 
mensions as well as in one dimension, due to the similar¬ 
ity of the main e quation s (see below) with the Gutzwiller 
approximatiorP^ 18 * 20 * 24 *, one can say the two approaches 
should succeed or fail together. Unfortunately, in two 
or three dimensions — which are cases of common in¬ 
terest — there is no a priori way to know whether the 
approximation will be good or poor. 

We note that when the decoupling approximation is 
performed and the constraint imposed only on average 
in Eq. (131, we may decide to change N m i n and N max to 
allow additional unphysical states with negative or posi¬ 
tive occupations. For example, letting N m i n —» — oo and 
N ma x —► +oo, which in turn makes Ci a irrelevant, yields 
the mean-field 0(2) slave-rotor formalism used in pre¬ 
vious worlP^DHl As noted previously®, allowing these 
unphysical states is not a major error in the limit of 
strong interactions since number fluctuations to unphys¬ 
ical occupancies are at any rate unlikely due to their 
large Coulombic penalties; however, for small interaction 
strengths, the unpliysical state have significant weight 
in the wave function which creates incorrect behavior, 
e.g., improper behavior of the quasiparticle weight ver¬ 
sus interaction strength. At the other extreme, a separate 
slave boson for each spin+orbital combination ima gives 
N m i n = 0 and N max = 1 which is just the “slave-spin” 
formalism ] 20 * 21 ! 
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With this separability assumption, the time- 
independent Schrodinger equation for the original 
system separates into two separate equations where 
the constraints of Eq. (13) are enforced by Lagrange 
multipliers hi a appearing in the two Hamiltonians. In 
the remainder of this section, we discuss the simplest 
U = U' and J = 0 case for simplicity. Full expressions 
involving U, U' and J for various slave boson choices 
follow after this section. The spinon Hamiltonian is 
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The spinons are coupled to the slave bosons via the aver¬ 
age {0\ a 0i' a ')s which renormalizes spinon hoppings be¬ 
tween sites i and i'. The spinon problem is one of non¬ 
interacting fermionic particles with spin. 

The slave boson Hamiltonian takes the form 
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where the spinon average ( fj mrJ fi'm'cr )/ renormalizes the 
slave boson hoppings. The slave boson problem is one of 
interacting charged bosons without spin. 

The original problem has been reduced to a set of 
paired problems that must be solved self-consistently. 
The spinon and slave boson problems only communicate 
(i.e., are coupled) via averages which renormalize each 
other’s hoppings. At this point, one must make some ap¬ 
proximations in order to solve the interacting bosonic 
problem. Typical approaches to date include single¬ 
site mean field approximationi^^, multiple-site mean 
fielcP^, approximation by sigma models to yield Gaus¬ 
sian integrals 17 ^as well as a combination of using tight- 
binding parameters obtained using Wannier functions 
from DFT followed by a mean-field approximation^. 

Separately, a procedure is needed to obtain the Ct a - 
To this end, at U = U' = J = 0, one chooses the Ct a 
to ensure that the spinon bands reproduce the original 
non-interacting band structure and associated occupan¬ 
cies (i.e., fillings). This means that the slave-boson ex¬ 
pectations {0\ a 0t’a 1 ) s should be unity in order not to 
modify the spinon hoppings away from the original non¬ 
interacting hoppings timi'm'a ■ The numbers Ci a and hi a 
are determined by making (0\ a 0i' a >) s unity as well as 
reproducing the non-interacting occupancies or Filings. 
This actually requires us to solve the coupled slave and 


spinon problems at U = U' = J = 0 self-consistently 
to obtain Ci a and hi a . The values of Ci a are then held 
fixed from that point forth when turning on U,U', J to 
non-zero values to self-consistently solve the actual inter¬ 
acting problem. 


Prior to solving some model problems within our new 
framework, we provide more complete descriptions of a 
number of potential choices for the slave-boson model 
(i.e., the choice of a) with full U, U' and J depen¬ 
dence. Differing choices split the interaction terms Hint 
of Eq. § in different ways between the spinon and slave 
sectors. This opens the door to systematic comparison 
between the different types of treatments of correlations 
with the slave bosons. 


1. Number slave 


The simplest approach is to simply create a single slave 
boson on each site i whose number operator Ni counts 
all the electrons on that site. In other words, the label a 
contains all the mcr orbitals on that site: it is superfluous 
so we can write Oi a = Oi. Description of the physically 
allowed states will set N m i n = 0 while N max will be the 
maximum number of electrons allowed on that site: e.g., 
10 for d shells or 14 for / shells. 


In this case, the slave boson can only represent the U 
term of the interaction in Eq. |2]); all remaining interac¬ 
tion terms must be treated at the mean-field level in the 
spinon sector. The slave Hamiltonian in this case is 
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while the spinon Hamiltonian contains all the remaining 
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interaction terms at mean-field level: 
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In the derivation of the expression for the above spinon 
Hamiltonian Hf , we have used the definition of the one- 
particle density matrix 

Pba = (fa fb) f > 

the standard mean-field contraction of four particle op¬ 
erators into two-particle operators weighed by averages 

fafbfcfd ~ Pdafbfc ~ Pcaffrfd + Pcbfafd ~ Pdbfafc > 

and the average occupations 


boson. We call this the orbital slave method. Here the 
index a labels a specific spatial orbital m and ranges over 
the two spin directions for that orbital: we have O im for 
the raising/lowering operator and N irn for the particle 
count slave operators. The slave sector can now directly 
describe more of the interaction terms: 
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and the spinon Hamiltonian is less complex than the pre¬ 
vious case as it only has the J terms at mean-field level: 
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3. Spin slave 


This approach has the simplest slave Hamiltonian and 
the most complex spinon Hamiltonian. This is because 
the number-only slave boson can only describe the sim¬ 
plest U part of the interaction; the remaining terms in¬ 
volving U' and J must be handled at mean-field level by 
the spinons. As mentioned above, the physical range for 
the occupation numbers of the number slave N-i is from 
zero to the physically allowed N max for that site. How¬ 
ever, we can decrease N min below zero and N max above 
the physical value if desired: in the limit where the range 
of occupancies allowed is very large, we automatically re¬ 
cover the 0(2) slave-rotor method. 

2. Orbital slave 

A more fine-grained model is to count the number of 
electrons in each spatial orbital m separately with a slave 


An alternative fine-graining is to have two slave bosons 
per site that count spin up and spin down electrons sep¬ 
arately but with no orbital differentiation. Namely, a 
labels a spin state a but ranges over all spatial orbitals. 
Hence, we have O ia and N la for our slave operators. The 
slave-boson Hamiltonian is 

*. = |e((e»..) 2 -eE 
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while the spinon Hamiltonian is 
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4- Spin+orbital slave 


This approach represents maximum fine-graining 
whereby we use a slave boson for each spin+orbital com¬ 
bination. Thus the index a now represents a full set of 
quantum numbers me r so we have Oi m(7 and Ni ma for the 
slave operators. The physically allowed occupancies are 0 
and 1 which is isomorphic to a pseudo-spin. For this rea¬ 
son, the name used for this approach in the literature is 
the “slave-spin” methocP^ii. However, given the possi¬ 
ble confusion this term creates between the real electron 
spin as well as the difficulty of using such a name un¬ 
ambiguously in our generalized formalism, we prefer the 
more explicit name “spin+orbital slave” where the spin 
refers to the physical electron spin. 

In this approach, we can describe the maximum num¬ 
ber of interaction terms in the slave Hamiltonian: 


U 


£« = oE E*M -EE 


U' — U +y 

+ E E^» 

m^m' \ cr 

j 


EE 


__/V- _/V- 

2 / j / j ly ima iy u 
a m^m' 

+EE hima-Nimo 


i ma 

ft 


^ ' timi'm'aifimcrfi'Tn'cr)fd ijner Oii m i a . (22) 
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spin flip and pair-hopping terms: 
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ima i m 


‘im'i'im 


El Pim'-'i'm' 


f Or'r 


' cr fimerfi' 


(23) 


We mention that in prior work 2 ^ the spin flip and 
pair hopping terms were argued to be well treated in 
the slave-particle sector instead. Namely, they were re¬ 
moved from the spinon Hamiltonian and the following 
terms were added to the spin+orbital slave Hamiltonian: 


^ y 1 

m^m' 

+ ^tm^tmi^im'^im'i + h ' C -) (24) 

where the S operators in the number basis are 

s+ =(;i!) • s_ =(sj)- < 25 > 

While such an ad hoc approach is not the strictly the¬ 
oretically consistent way to split operators between the 
spinon and slave boson sectors, in practice it does re¬ 
produce the desired behavior of the spin flip and pair 
hopping terms in the slave boson sector and does not 
introduce any numerical difficulties. 

Our approach provides some insights into this incon¬ 
sistency issue while simultaneously easing some techni¬ 
cal problems that can arise in the spin+orbital slave ap¬ 
proach. Part of the inconsistency is that the ad hoc S ± 
slave operators are not the same as the O operators (the 
only way to make them the same is the extreme choice 
Cimcr = 0). For example, at half-filling when Ci ma = 1, 
6 and Of are the same operator and equal the S x Pauli 
matrix, so one can not use O and O' to represent any sen¬ 
sible representation of the the spin-flip or pair-hopping 
interaction terms. If we choose to include some unphys¬ 
ical states, however, things become different. For exam¬ 
ple, if we widen the range of the spin+orbital slave bo¬ 
son occupancies from {0,1} to { — 1,0,1,2} then O and 
O become different. One can then write a more natural 
interaction term of the form 

— J E ^\m4^imiO\ m ,^O im i^ 

m^m' 

+ ^imt^imi^ im 't^im'i + h.C.) (26) 


The corresponding spinon Hamiltonian still contains the 






that only uses the slave boson O operators. 

Separately, enlarging the set of occupancies beyond 
{0,1} has some technical advantages. When the occu¬ 
pancies are limited to {0,1}, then we have the analytical 
forirPO given by C ima = [n im(T { 1 - n irncr )\~ 1/2 . For oc¬ 
cupancies nearing the extremes (riima -+ 0 or 1), the 
Cima become very large, the Oim a matrices become ill- 
behaved, and the numerical algorithm using them be¬ 
comes difficult to stabilize. Permitting a wider range of 
occupancies such as { — 1,0,1,2} makes the Ci ma have 
reasonable values and the numerical procedure is well 
behaved. In the large interaction limit, the addition of 
the unphysical states (e.g., occupancies —1 and 2 in this 
case) is not a major error since number fluctuations are 
suppressed; the main problem is in the weak interaction 
regime where this enlargement of the occupancy basis 
provides numerical stability but can produce quantita¬ 
tive errors. 


III. MEAN-FIELD TESTS 

We now proceed to describe computational results 
based on a simple single-site, paramagnetic, nearest- 
neighbor, mean-field solution of the slave Hamiltonian at 
half filling. This will permit us to both reproduce prior 
literature as well as to compare various slave Hamiltoni¬ 
ans to each other. 

To do so, we shift the local interaction energies so that 
they are zero at half filling, i.e. when ro imtT = 1/2. We 
also make the standard choice U' = U — 2J. The local in¬ 
teraction term (ignoring for the moment the spin flip and 
pair-hopping terms) takes the form from prior work 20 : 

TJ _ O T 

Hint = -o-^ “ nl ° rh)2 + J E^ im ~ X ) 2 

(27) 

- ^ ~ n orb / 2 ) 2 

a 

where n l orb is the number of localized correlated spatial 
orbitals on site i. 

In the single-site mean-field approximation, we will be 
solving for a single site self-consistently coupled to an 
averaged bath of bosons on the nearest neighbor sites. 
Our assumptions ensure that all sites are identical with 
no spin polarization. Furthermore, to connect to the lit¬ 
erature, we further assume that in the multi-orbital case 
there are only non-zero hoppings between nearest neigh¬ 
bor orbitals with the same to index. With all these as¬ 
sumptions, it is easy to see that C/ a = 1 is the choice that 
gives half-filling for the slave problem at U = U' = J = 0. 
In addition, we can set the Lagrange multipliers hi a = 0 
since we have set the half-filling energy to be zero. The 
density matrix elements ( fimafi'm’a)f that renormalize 
the slave boson hoppings will be spin and site indepen¬ 
dent and will be non-zero only when to = to’. Hence, 
they can be absorbed into the definition of the hopping 


elements ti m i> m t a - The density of states for a spinon band 
is taken to be the standard semicircular one 


D(E) 


,V4t 2 -E 2 


2nt 2 


o 


if |f;|< 2\t\ 
else 


We begin with J = 0. The slave Hamiltonian is 

Hs — ^ ^ l Hi a Tl 0 rb 

i \ a 

"EE (28) 

iot ra£a 

where the effective hoping for spatial orbital to. is 

^ = E E timi'm'olPa') s . (29) 

i' ol' m' g^cx' 

The simple form of this Hamiltonian makes it easy to 
directly read off the quasiparticle weight renormalization 
Z a which narrows the spinon bands: 

= (O a ) 2 s . (30) 

When Z a = 0, a Mott insulator is realized in such a 
simple single-site model 27 ! We solve the problem self- 
consistently for different slave models. Since at half¬ 
filling the Lagrange multipliers /i,; a = 0, all that is re¬ 
quired to solve the spinon problem is to renormalize each 
spinon band width (i.e., hopping) by the appropriate Z a 
factor. 


A. Single-band Mott transition 

We begin with a single-band model where there is one 
spatial orbital per site. Figure [I] compares various slave 
models based on the dependence of Z on U. Specifically, 
we compare the slave rotor model (allowed occupancies 
from —oo to +oo), the orbital slave model (allowed occu¬ 
pancies 0, 1, or 2) which here is identical to the number 
slave model, the spin+orbital slave (“slave-spin”) model 
(allowed occupancies 0 or 1) and the Gutzwiller approx¬ 
imation where Z Gu tzwiUer = 1 - ( U/U c ) 2 . 

For this system, the Gutzwiller and spin+orbital slave 
methods predict exactly the same results, as noted 
previouslypSl In fact, the spin+orbital slave model, at 
half-filling for a single orbital per site at the single¬ 
site mean field level, can be shown to be isomorphic to 
the Gutzwiller approximation as well as to the Kotliar- 
Ruckenstein model as described by Bunemann.^ This 
shows that, beyond their utility as mathematical mod¬ 
els, such slave-boson methods can parallel and help un¬ 
derstand other approaches that originate from apparently 
different sets of many-body approximations. 

The slave-rotor method has an aberrant behavior for 
small U. Specifically, Z for the slave-rotor method has 
the small U expansion 

Zrot.or = 1 — 0{\ju/t e ff) . 


(31) 
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The reason for this behavior is due to the unbounded 
number states permitted in the slave rotor model. Specif¬ 
ically, in the number basis the slave-rotor problem corre¬ 
sponds to an infinite one dimensional lattice labeled by 
Ni, with hoppings t e ^ between neighboring sites, and 
with a quadratic potential UN?/ 2. For small U, the 
ground state of this problem will be spread over many 
sites so that we can take the continuum limit. The prob¬ 
lem turns into the textbook one dimensional harmonic 
oscillator with mass 1/(2 t e ff) and spring constant U. 
The ground state wave function ip(Ni) is a Gaussian, and 
( 0) s = '}Z n ’4 , {n)ip{ri — 1) can be computed. Expansion 
in U then gives the above form. 

In reality, however, perturbation theory guarantees 
that quasiparticle weights are modified starting at sec¬ 
ond order in the interaction strength: 

Z = 1 - 0((U/t eff ) 2 ) . (32) 

The slave-rotor fails since for small U it spreads the wave 
function over a large number of unphysical states. What 
this means is that one would incorrectly overestimate 
the importance of electronic correlations at weak inter¬ 
action strengths when using the slave-rotor method. In 
this view, our orbital and number slave methods may 
be viewed as corrected rotors which are restricted to the 
appropriate finite set of physical states. Finally, Fig- 
ure[l]illustrates that slave methods employing finite slave 
Hilbert spaces all automatically correct the small U be¬ 
havior. 


B. Isotropic two-band Mott transition 

Next, we consider a two-band degenerate Hubbard 
model. Figure [2] displays the results. We note that the 
two band e g model is of physical relevance as the slave- 
rotor has shown itself to be of use in e g nickelate systems 
within a pd modeP-^. For this particular degenerate case 
with high symmetry, the spin slave and orbital slave mod¬ 
els turn out to be identical since each posits two slave 
particles each with the allowed occupations 0, 1, or 2. 
We note that, in this case, the slave rotor and number 
slave become very similar for large U: once slave number 
fluctuations of N t are small, the size of the slave Hilbert 
space becomes irrelevant. 


C. Anisotropic Orbital-Selective Mott Transition 

We present mean-field calculations exemplifying the 
orbital-selective Mott transition in an anisotropic two 
band model with paramagnetic solution and at half fill¬ 
ing. We take spatial orbital m = 1 to have the larger hop¬ 
ping ti while m = 2 has the smaller hopping t 2 - Hence, 
h/ti specifies the degree of anisotropy. 

The first slave model for this system is the s pin+orbital 
method which has been used previouslj® 2 fl each slave 



FIG. 1. Quasiparticle weight Z as a function of U/U c for 
different slave-particle models for the paramagnetic single¬ 
band Hubbard at half filling. U c is the critical value of U when 
Z = 0, i.e., the Mott transition, for each model. The black 
squares show slave rotor results, the blue line is the Gutzwiller 
approximation results which for this model are the same as 
the spin+orbital slave (“slave-spin”) results in blue crosses, 
and the red circles show the orbital slave results (identical 
to the number slave). We note that the orbital-slave Hilbert 
space is very small, so that it does not agree with the rotor, 
unlike the two-band slave number. 



FIG. 2. Quasiparticle weight Z as a function of U/U c for 
different slave-particle models for a degenerate paramagnetic 
two-band Hubbard model at half filling. 


boson has allowed occupancies 0 or 1. The second model 
is to forgo the explicit spin degree of freedom in the slave 
description and to employ the orbital slave model where 
each slave boson has allowed occupancies 0, 1, and 2. 
The comparison tests the importance of explicit treat- 
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ment of spin in the electronic correlations for such a sys¬ 
tem. We will focus on the Orbital-Selective Mott Tran¬ 
sition (OSMT) when one orbital has a finite bandwidth 
and is metallic while the other has undergone a Mott 
insulating transition and is localized. 

We begin with J = 0. Figure [3] illustrates the behavior 
of the renormalization factor Z for both bands versus 
U for three different t 2 /t\ ratios within the two slave 
particle models. An OSMT occur for small enough t 2 /t\ 
ratio but the critical value depends on the type of slave 
model. For the orbital slave model, we find that OSMT 
occurs when t 2 /t\ < 0.25 while for spin+orbital slave we 
must have a slightly smaller value of t 2 /t\ < 0.2. 

We now consider J > 0. We continue to treat the 
spinon problem as that of a simple, paramagnetic, half- 
filled tight-binding model with two separate bands with 
each hopping renormalized by the appropriate ( O a ) s . For 
the orbital slave model, we can only include the first two 
terms of Eq. (27) due to the lack of an explicit spin label 


in the slave description. Thus we will compare the orbital 
slave and spin+orbital slave using the same interaction 
term 

Hi nt = U ^(N i -2) 2 + jY / (Nim-l) 2 . (33) 


It is clear from the above two interaction terms that, for 
fixed C/, J > 0 permits larger orbital independent num¬ 
ber fluctuations (i.e., it reduces the correlation effect of 
this mode) since U' = I/ — 2J becomes smaller in the first 
term. However, the second + J term simultaneously pun¬ 
ishes intra-orbital number fluctuations and thus enhances 
intra-orbital correlation effects which in turn favors an 
OSMT. 

The phase diagram as a function of t 2 jt\ and J for 
this system in shown in Figure [4j The boundaries shown 
separate regions where OSMT occurs (above the bound¬ 
aries) from where a standard Mott transition occurs (be¬ 
low the boundaries). The figure confirms the fact that in¬ 
creasing J favors OSMT. Qualitatively, the orbital slave 
and spin+orbital slave show very similar behavior: they 
have a critical t 2 /ti at J = 0 between 0.2 — 0.25 for 
OSMT to occur, and then with increasing J the criti¬ 
cal t 2 /t\ becomes larger so less anisotropy is needed to 
drive an OSMT, as observed previously in DMFTp3 and 
spin+orbital slave calculation^®!. 

We have also considered the case where we permit 
the orbital slave model to have unlimited occupations: 
namely, we have a two rotor model (one for each orbital 
occupation). In this case, we find that no OSMT is possi¬ 
ble when J = 0 for any bandw idth ratio t 2 /t\. This result 
is similar to previous DMFftf^SlM^ which found that a fi¬ 
nite J is needed in order to have an OSMT. However, it 
disagrees with the results of previous orbital+spin slave 
results®! as well as our results above where we find that 
a small enough bandwidth ratio t 2 /t\ makes an OSMT 
possible even for J = 0. These differences further il¬ 
lustrate the need for multiple models and cross verifica¬ 
tion when describing a possible OSMT for real materials 



FIG. 3. Quasiparticle weights for the paramagnetic 
anisotropic two-band single-site Hubbard model at half filling 
as predicted by the orbital+spin slave model (blue squares) 
and the orbital slave model (red circles) at J = 0 for three 
t 2 /t\ ratios. In each plot, the Z value for the first orbital with 
larger hopping t\ is denoted by symbols and a continuous line 
while for the second orbital symbols and a dashed line are 
used. An OSMT occurs when the two Z do not go to zero at 
the same U value: orbital slave (red circles) in the center plot 
and both slave models in the lower plot. 


which have complex band structures (e.g., the three-band 
Ca 2 _a;Sr a ,Ru 04 systenP). 

Prior woi'k^has shown that the presence of the Hund’s 
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term 


„ ^ ' 'y ' (Pma l/2)(n TO ' cr 1/2) 


(7 


= -^ EI (^ -1 ) 2 ■ ( 34 ) 


makes OSMT slightly more difficult to achieve as it in¬ 
creases inter-orbital m ^ m 1 correlations by favoring spin 
pairing between different orbitals but does not aid intra¬ 
orbital correlations. Separately, adding the spin-flip and 
pair-hopping terms makes OSMT easier to achieve^. 

Although not directly relevant to our main focus, for 
completeness we include a final comparison based on a 
fixed slave model with various combination of interac¬ 
tion terms. We choose the spin+orbital orbital model 
and then choose to include different interaction terms in 
the slave-particle Hamiltonian. The first choice is the 
interaction terms used above in Eq. (33). The second 
choice is to add the Hund’s term: 


TT _ O J 

Hint = - < rb ) 2 + - l ) 2 




(35) 


Prior work 2 ^ has shown that the presence of the Hund’s 
term makes OSMT more difficult to achieve as it in¬ 
creases inter-orbital correlations by favoring spin-pairing 
among different orbitals but does not enhance intra¬ 
orbital correlations. 

The third choice is to add the spin-flip and pair¬ 
hopping terms as per the ad hoc method of Eq. (24): 


JLTl _ 

**int 


U-2J 


a 

^ y y ( E m t E m' v E m '" 

m+m' 

^ irojE im' 4- h.C.’) . 

Adding these spin-flip _and pair-hopping terms makes 
OSMT easier to achievd^. 

Phase diagrams for the second and third choices above 
are available in the literature 20 and are reproduced in 
Figure [5] which also includes the results of the first choice 
as well. We note that only including the intra-orbital 
terms (first choice) or all terms (third choice) leads to 
essentially the same phase diagram. However, exclud¬ 
ing the spin-flip and pair-hopping terms (second choice) 
makes it harder to achieve an OSMT phase: one can not 
achieve an OSMT for any reasonable J once the band¬ 
width ratio ti/ti exceeds ~ 0.6. The physics behind this 
progression is as follows. Starting with J = 0 and a rel¬ 
atively large U, the ground-state basically contains only 


(iV i -2) 2 + 


(36) 



FIG. 4. Phase diagram for the anisotropic two-band single¬ 
site Hubbard model at half-filling as a function of the 
anisotropy ratio O/ti and J. Two slave boson methods 
are used: orbital slave (red circles) and spin+orbital slave 
(blue squares). In each case, the boundary curve demarcates 
the possible existence of an Orbital-Selective Mott Transition 
when U is ramped up from U = 0. Regions above the bound¬ 
ary display OSMT while regions below it present a standard 
Mott transition where both bands become insulating at the 
same critical U c value. 


states which are half-filled and have a total of two elec¬ 
trons per site (there are six such states). Adding the 
intra-orbital term (first choice) with J > 0 then further 
restricts us to the four states with only one electron per 
orbital (but with no preference for spin states). Such a 
ground-state can suffer an OSMT when further increas¬ 
ing U since the narrower band (more localized orbital) 
can become fully localized. Next, adding the Hund’s 
term (second choice) creates a preference for the two 
spin-aligned states in this four dimensional subspace by 
lowering their energy: this enhances inter-orbital corre¬ 
lations at the expense of intra-orbital correlations which 
favor an OSMT phase. Third, adding the spin-flip and 
pair-hopping (third choice) terms essentially cancels the 
effect of the Hund’s term. This is explained by a straight¬ 
forward computation of the matrix elements of this inter¬ 
action in the four dimensional subspace. One finds that 
the spin-flip term couples the two states where electrons 
have opposite spins with a strength that is precisely such 
that their symmetric combination has the same energy 
lowering as the Hund’s term induces for the spin-aligned 
states. Thus, we are essentially back to the four states 
we had when only operating with the intra-orbital in¬ 
teraction (first choice). Our final comment is that these 
differences are not very dramatic once the hopping ratio 
t 2 /t\ is below ss 0.5. As Fig. [5] shows, in all cases only 
a modest value for J is sufficient to stabilize the OSMT 
phase instead of a standard Mott transition. 
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FIG. 5. Phase diagram for the anisotropic two-band single¬ 
site Hubbard model at half-filling as a function of the 
anisotropy ratio tilt\ and J for the spin+orbital slave model. 
Three different interaction terms are used: intra-orbital term 
only which is Eq. ( |33| ), intra-orbital plus Hund’s which is 
Eq. (351, and all terms included which is Eq. (361. 


D. Ground State Energies 

A final and most stringent test for the slave models 
is to compare their total energies. In the interest of 
space, we will focus on the simplest case of degenerate 
orbitals, isotropic hopping, and phases that are para¬ 
magnetic and paraorbital (no orbital differentiation) to 
make some general comments. In a fully self-consistent 
model with more parameters and non-degenerate bands, 
we may expect more complexity to be revealed. Previous 
workSD has shown that ground-state calculations can re¬ 
veal competition between the orbital-selective Mott state 
(due to very large crystal-field splitting) and an anti¬ 
ferromagnetic Mott insulating state (due to a large J), a 
transition which is likely first-orderli. 

With J = 0, the ground state energy per site of the 
paramagnetic and paraorbital phase is 

Eg = ~ E E - n °rtf) ■ ( 37 ) 

a ct a 


the high degree of spin and orbital symmetry since the 
Hartree-Fock Slater determinant wave function will be 
unchanged versus U and always predicts a metallic sys¬ 
tem. 

The next observation is that for small U, some of the 
slave models do worse than Hartree-Fock. However, as U 
is increased their total energies eventually drop below the 
Hartree-Fock one. Furthermore, increasing the number 
of bands from one to two improves the total energies of 
all slave methods compared to Hartree-Fock. For a given 
number of bands, increasing the fine-grained of the slave 
model (i.e., having more slave modes per site) also lowers 
the total energy. Hence, the slave-rotor is generally the 
worst performer. 

A final observation is that only the fully fine-grained 
spin+orbital slave method always predicts a total energy 
below that of Hartree-Fock. It also has the correct linear 
slope of E g versus U matching the Hartree-Fock one. The 
other slave methods have higher slopes of E g versus U at 
the origin so that they can only outperform Hartree-Fock 
beyond some finite value of U. The slope matching of the 
spin+orbital slave is a natural expression of its account¬ 
ing in detail for all the quantum numbers on each site 
and in being forced (like all slave models) to reproduce 
the non-interacting state at U = 0. The fact that the 
other slave models have higher slopes is a reflection of 
their larger (and quantitatively incorrect) number fluc¬ 
tuations at U = 0. Namely, the interaction Hamiltonian 
H mt is a quadratic function of the occupancy numbers 
so that its expectation value (the interaction energy) de¬ 
pends directly on the fluctuations of these occupancies; 
at fixed U , the larger the set of allowed occupancies in a 
slave model, the larger this quadratic fluctuation and the 
higher the interaction energy. In fact, the number fluctu¬ 
ations of the slave-rotor model are so large at U = 0 that 
they lead to a pathological infinite slope of E g versus U 
at U = 0. By comparison, the number slave method, 
which can be viewed as a corrected rotor, has a much 
more reasonable behavior. 

As a side note, it is interesting that for the single¬ 
band case, one has the following analytical results based 
on the coincidence of the of the spin+orbital slave and 
Gutzwiller approximations. In the metallic phase, where 
U < U c , the quasiparticle weight Z is given by 

Z = 1 - U 2 /U 2 (38) 


We compute the ground-state energy as a function of U 
for one-band and two-band isotropic models at half-filling 
(same systems that are in the above sections) and also 
include the Hartree-Fock total energy. Figures [6] and [7] 
display the energies versus U for the one-band and two- 
band cases, respectively. The plots employ the half-band 
width D = 2 1. 

In all cases, for large enough U the slave models pro¬ 
duce an insulating phase (i.e., isolated atomic-like sites) 
which has zero hopping and zero number fluctuation and 
thus zero energy in this model. The Hartree-Fock to¬ 
tal energy necessarily has a linear dependence on U for 


and from perturbation theory at small Z 18 
U C = 8D. 

Using the following definition: 

to = t{f ima f ima ) U=0 
the ground-state energy is given by 


U 


U 2 


Eq 2to + o . 

9 ° 4 128tn 


(39) 


(40) 


( 41 ) 
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FIG. 6. Ground-state energy per site ( E g /t ) of a single band 
Hubbard model at J = 0 in the paramagnetic phase at half 
filling for a variety of slave representations as well as for the 
Hartree-Fock approximation. D = 2t is the band width of 
the non-interacting system. For this model the orbital slave 
is identical to the number slave and the spin slave is the same 
as the spin+orbital slave. 
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FIG. 7. Ground-state energy per site ( E g ) for an isotropic 
two-band Hubbard model at half filling for J = 0 in the para¬ 
magnetic and paraorbital phase. 


fine-grained model which is in the spin+orbital slave rep¬ 
resentation (“slave-spin” in the literature). 


IV. CONCLUSIONS 

We have developed a generalized formalism that repro¬ 
duces previous slave-particle formalisms in appropriate 
limits but also allows us to define and explore interme¬ 
diate models and to compare them systematically. Our 
formalism moves beyond the analogy with angular mo¬ 
mentum behind slave-rotor formalism, and instead works 
directly in the physically correct finite-sized number rep¬ 
resentation permitting new models to be developed in a 
more natural way. As an example, we have shown how 
the standard Mott transition as well as the orbital selec¬ 
tive Mott transition appear in different slave models for 
single-band and two-band Hubbard models. 

We believe it is useful to have a variety of slave par¬ 
ticle methods on hand as they provide computationally 
inexpensive methods for exploring the role of electronic 
correlations in materials and interfaces with broken sym¬ 
metries (e.g., orbital symmetry breaking). The cheap 
computational load is particularly advantageous for in¬ 
terfacial systems where translational symmetry is lost 
in one direction and simulation cells that capture the 
region near the interface must contain at least tens to 
hundreds of atoms. As such, these simpler slave-particle 
models are useful for exploratory research where more 
accurate and expensive Hubbard-model solvers such as 
DMFffEMl 

would be prohibitive to apply routinely. The 
ability to isolate potentially interesting correlated degrees 
of freedom from each other by choosing different slave ap¬ 
proaches may illuminate which degrees of freedom are the 
most critical to model accurately. 
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For the insulating state (U > U c ), we have E g = 0. 

Our calculations in this section permit us to say that 
while our generalized approach permit us to easily com¬ 
pare different slave models and isolate different degrees 
of freedom simply, total energy comparisons are much 
more challenging. First, one should do energy compar¬ 
isons of different phases within a single slave model since 
the differing models can produce differing total energies 
with dependence on the details of the system. Second, 
after understanding the relevant degrees of freedom and 
how the influence the physical behavior, the total en¬ 
ergy calculation should be most accurate with the most 


Appendix A 

In this appendix, we provide some detailed examples of 
how the physical subspace is isolated from the extended 
Hilbert space of spinon+slave boson states and how the 
operators act in the physical subspace. In the process, 
we also provide explicit examples for various choices of 
the slave labels a. We focus on a single site i and hence 
suppress the site label i below. 

The original Hilbert space, i.e, the Fock space of the 
fermionic d mcr field operators, is spanned by basis kets in 
the occupancy number representation for the field opera¬ 
tors and have the form |{n m(T }) where n ma £ {0,1}. The 
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enlarged Hilbert space for spinons and slave particles is 
spanned by product kets in the number occupancy basis 
of the form 


An orbital slave representation has slave bosons counting 
the number of electrons in each spatial orbital, so a = m. 
The 16 allowed kets are 


|{foncr})/|{-^a}) s 

where, again, n ma € {0,1} are the fermionic spinon oc¬ 
cupancies while N a are the bosonic particle counts. The 
/ and s subscripts label the spinon and slave boson kets. 

The constraint of Eq. ^ on the physical allowed states 
translates to the numerical constraint 

N a = E ^ma • (Al) 

ma^ct 


10,0,0,0^10,0),, 11, 0,0,0) / |l, 0) s , 10,1,0,0^11,0),, 
10,0,1,0^10,1),, 10,0,0,1^10,1),, 11,1,0,0^12,0),, 
11,0,1,0^11,1),, 11,0,0,1)/| 1,1), , 10,1,1,0^11,1),, 
|0,1,0,1)/11,1), , 10,0,1,1^10,2),, 11,1,1,0^12,1),, 
11,1,0,1^12,1),, 11,0,1,1^11,2),, 10,1,1,1^11,2),, 

11 , 1 , 1 , 1 ^ 12 , 2 ),. 


We remember that we choose the {n ma } to match exactly 
between the original electron and spinon kets. 

We begin with the simplest example of a single 
spatial orbital on the site where the kets look like 
\n^,ni) f\{N a }} s . There are two states for electrons and 
thus a total of four possible configurations: no electrons, 
one spin up electron, one spin down electron, and a pair 
of spin up and down electrons. If we have a single slave 
boson per site to simply count the number of electrons so 
the a label is nil (i.e., the number slave representation), 
then our four physically allowed kets are 

|0,0) / |0>, , | 1 , 0 ) / | 1 >, , |0,1> / |1>, , 11 , l)j|2), . 

We note that the number of slave particles is constrained 
by Eq. © to the total the number of spinons. 

Next, if we have this single orbital but instead we 
choose to have a slave mode per spin channel (i.e, the 
spin+orbital slave representation), then we have two sets 
of slave bosons since now a = a. The four physical states 
are now 

| 0 , 0 ) / | 0 , 0 >, , | 1 , 0 > / | 1 , 0 ), , 10 , 1 ^ 10 , 1 ),, 11 , 1 )^ 11 , 1 ),. 

A more complex set of examples has two spatial or¬ 
bitals per site. Here we have four choices of spinon label 
ma which we order as If, 1|, 2f, 24,. For the number 
slave representation, we have the 16 physical kets 

10 , 0 , 0 , 0 ^ 10 ), , | 1 , 0 , 0 , 0 > / | 1 >, , 10 , 1 , 0 , 0 ^ 11 ), , 
10 , 0 , 1 , 0 ^ 11 ), , 10 , 0 , 0 , 1 )^ 11 ), , 11 , 1 , 0 , 0 ^ 12 ), , 
11 , 0 , 1 , 0 ^ 12 ), , 11 , 0 , 0 , 1 )/ 12 ), , 10 , 1 , 1 , 0 ^ 12 ), , 
10,1,0,1^12), , 10,0,1,1)^12), , 11,1,1,0)y|3), , 

11,1,0,1)/|3), , 11,0,1,1)/13), , 10,1,1,1^13), , 

11 , 1 , 1 , 1 )/| 4 ), • 


Alternatively, one can use the spin slave representation 
where the bosons count the number of electrons of each 
spin, so a = a. The allowed kets are 


10 , 0 , 0 , 0 ^ 10 , 0 ),, 11 , 0 , 0 , 0 ^ 11 , 0 ),, 10 , 1 , 0 , 0 )^ 10 , 1 ),, 
10 , 0 , 1 , 0 ^ 11 , 0 ),, 10 , 0 , 0 , 1 ^ 10 , 1 ),, 11 , 1 , 0 , 0 ^ 11 , 1 ),, 

11 , 0 , 1 , 0 ^ 12 , 0 ),, 11 , 0 , 0 , 1 )^ 11 , 1 ), , 10 , 1 , 1 , 0 ^ 11 , 1 ),, 

10 , 1 , 0 , 1 ^ 10 , 2 ),, 10 , 0 , 1 , 1 )^ 11 , 1 ), , 11 , 1 , 1 , 0 ^ 12 , 1 ),, 

11 , 1 , 0 , 1 ^ 11 , 2 ),, 11 , 0 , 1 , 1 ^ 12 , 1 ),, 10 , 1 , 1 , 1 ^ 11 , 2 ),, 

11 , 1 , 1 , 1 )^ 12 , 2 ),. 


The final point is to check that the original electron op¬ 
erators d ma have the same effect as the combination of 
spinon and slave f m aO a in the physical subspace. That 
this is in fact true follows directly from the defining Equa¬ 
tions (7|9) along with the constraint on N a in Eq. (All. 


It is easy to check that the matrix elements of d ma and 
fmaO a must match: 


{ n ma\dmcr\ n mcr) — / ( n ma \fmtr l^mcr) f ' s(N a \O a \N a ) g . 


The matching of the d and / matrix elements is clear be¬ 
cause the occupancies n mcr and n' ma match by definition 
on both sides and both operators have identical behav¬ 
ior on the occupancies as per Eqs. 0 and (|8]). Thus 
both sides are non-zero only if the n' occupancy set has 
one fewer total count than the n occupancy set. As long 
as N a > 0, the matrix element of O a is unity because 
N' a = N a — 1 m ust be true due to the occupancy match¬ 
ing of Eq. (Al). If N a = 0, it must be that n ma = 0, so 
that the matrix element of O a is irrelevant because the 
fermionic matrix elements (of d and /) are both zero. 
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